Oscillator glass in the generalized Kuramoto model: 
synchronous disorder and two-step relaxation 
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We consider a very general form of the Kuramoto model (KM) and derive equations for the 
macroscopic parameters of its stationary states. Surprisingly, we discover that a class of simple, 
analytically treatable KM variants can exhibit the long-sought oscillator glass behavior, which rep- 
resents a kind of synchronous disorder. We derive exact conditions for the appearance of the glassy 
states and study their exotic properties. Our findings offer the possibility of creating such states in 
real systems described by the KM, which may lead to a diversity of interesting phenomena. 



The Kuramoto model (KM) [l| was developed to pro- 
vide an analytically tractable description of the large 
populations of coupled oscillators that appear in many 
branches of science (HQ. It was used e.g. to describe the 
collective behavior of lasers 0, 0] , neurons in the brain 
@, Josephson junction circuits Q, even humans Q. We 
consider the KM in its very general form 



k- N 



ft -7j)> 



(1) 



where 8i(t) and w,; are respectively the ith oscillator's 
phase and natural frequency, N is their number, fc^ 
(A + 7j) represent the coupling strengths (phase lags) 
between the ith and jth oscillators, and parameters 
I\ = {u>i,ki,qi, /3i,7i} are drawn from a joint probabil- 
ity density G(T) = G(lu, k, q, (3,j). The case considered 
earlier 9] corresponds to (Q]) with q L = l,/3; = 7.; = 0. 

The KM has not previously been treated in such a gen- 
eral form ([TJ , although it includes the many KM modifi- 
cations and extensions studied earlier as particular cases. 
Thus, not much is known about possible behavior that it 
may exhibit. Daido reported evidence that the KM 
with random kiq^ and + 7, can undergo a glass transi- 
tion but, due to the complexity of the model, the resul- 
tant "oscillator glass" and its properties still remained 
mysterious I10l - ll4j . To illuminate this issue, some au- 
thors (e.g. [la]) attempted to find the same state in sim- 
pler models, but were not successful in this respect. 

In this Letter, we report our discovery of the oscillator 
glass in a class of simple, analytically treatable, models 
included within (JXJ) . We first generalize the framework of 
Q to encompass ([1]), and then proceed to consideration of 
the glassy states, deriving conditions for their appearance 
and studying their properties. 

Background. The oscillators' collective behavior in (p} 
can be described by two complex parameters 



N N 

Z = Re>* = -Y e"' , Y = We** = ±- V 



(2) 

where Z is the mean field whose amplitude R quantifies 
the degree of agreement between the oscillators' phases 
9i, while Y represents the weighted mean field; when 



qi = l,7j = 0, the fields coincide. The KM ([T]) can also 
be rewritten as 9i = uji — kiW J2j sin(#i — 9j — — $). 

In the continuum limit N — > 00, we drop indices and 
treat (JTJ) using the probability density function (PDF) 
f(9,T,t), representing the probability that the oscillator 
has parameters T and phase 9 at time t. The PDF can 
usually [IB-Gil be represented by the OA ansatz [2pj 



f(9,u,K,t) = 



G(T) 



2tt 



1 + 2Rc- 



a(r, t)e l 



(3) 



Following the conventional procedure [2 
that in our case a should satisfy 



! - a{T,t)e ie . 

one can show 



da 
~dt 



-<<"+*>) =0, (4) 
Z = Re 1 * = J a*(T,t)G(T)dT, (5) 



Y = We 



i<i> 



qe^a*(T,t)G(T)dT, 



(6) 



where dT = dkdqdfid"/ and integration is taken over 
the whole domain if unspecified [(—00,00) for k,q and 
[— 7r,7r] for (8,7]. Note that the difference $ — ^ repre- 
sents another meaningful parameter: it is invariant under 
6 ->• 8 — 9o (t), unlike $ and * individually. 

Parameter normalization. The model fTJ is invari- 
ant under both fcj — > — — > + tt and rescalings 
h — » ki/a,qi — > aq\. To avoid ambiguity, one can fix 
the normalization of % and specify a rule for choosing 
the sign of fcjj but different normalizations might be pre- 
ferred in different situations [2l|. Because the normal- 
ization can be fixed at any time, we will leave this ambi- 
guity to preserve generality. Note, that W (but not R) 
changes under the (fe, g)-rescaling and thus can be higher 
than unity, but that \k\W is invariant. 

Terminolog y a nd notation. We adopt similar termi- 
nology to @, |22|. The KM form ((TJ) is invariant un- 
der 9 = 8 — fit, which just changes the distribution 
G(ui, k, q, j3, 7) — s- G(ui + fl, k, q, /3, 7), with u> always de- 
noting the natural frequency in the current frame (with 
phases 9). Thus, one can consider the KM in different 
frames rotating at frequency f2 with respect to some fixed 
frame. For the latter we select the frame with zero mean 
frequency (u) = J ujG(T)dT = and call it the natural 



2 



frame; the distribution G(T) is denned in this frame. A 
stationary state (SS) is a state with a time- independent 
PDF |f = 0, which also implies Z = Y = 0. It can 
be stationary only in a particular rotating frame, so it is 
characterized by its frame frequency f2 and mean fields 
Z, Y. We call SSs with il = as natural states (NS), 
while SSs with O =/= are traveling wave (TW) states. 
For later convenience we define 

G^(T) = G(±uj + Q,k,q,p n ), 

Ti a \ A A p/a N ^VTT^/27r (7) 



i/ 2 +sin^(>/2) 



Note that P(</>,0) = 8{cf>) and P(</),oo) = l/2w. 
Additionally, we introduce 

I(u,k) = J e lf} G(T)dT, I± = J(±w + fi), 

j(u,k) = /" ge j(,3+7) G(r)dr, j± =/(±w + n). 



(8) 



As will be seen, I and J represent some kind of an "effec- 
tive" complex distribution of lj, k for the determination 
of Z and Y, respectively. 

Stationary states. We will treat SSs in frames where 
they are stationary, so that the parameter distribution 
becomes G(T) -> G^(T). By definition and ©, SSs sat- 
isfy ^2=0. Substituting this into (j4} and taking account 
of the OA validity condition |a| < 1 [20j . one finds that 
all SSs in their own rotating frames are described by 



a s e 



l($+/3) _ J V* a WCT-« , if H < \k\W, 



■ u — sign(aj) V^ 2 — k 2 W 2 



kW 



if 



> \k\W. 



(9) 

The other solution corresponds to an unstable position 
on the phase circle, and so is never realized. 

Using ([9|) in ((HJ , one can recover the PDF of an SS 



fs(e,r,t) 
G£(r) 



Vui 2 -k 2 W 2 /2TT 

|u-fciy sin(0-/3-*)| 



1 - - arcsin(^) - $ + nK(-k)) 
if \u\ < \k\W, 
if \u\ > \k\W, 

( 10 ) 

with H(-) denoting a Hcaviside function. 

Self- consistency and stability conditions. Using 
the self-consistency conditions (SCCs) onF ff become 



F 



dk 
kW 



\k\w 



vmv 2 



\k\W 



to 2 H^dT + i I x 



f Voj 2 -k 2 W 2 [J+ - J n }dT} = W. 

J\k\W " ' 



■tf dr- 



Taking its real and imaginary parts yields two equations, 
from which the SS parameters W, Q can be determined. 
They can then be used in (|9]),([5|) to determine Z as well. 

Next, applying the standard procedure [Hj], i.e. per- 
forming a linear stability analysis of ([3]) above incoher- 
ence (a = 0) and using a self-consistency argument (see 



@ and references therein), one can show that incoherence 
loses stability when there exists a solution Q to 



Trim / kJ(n,k)dk 
TrRe / kJ(tt,k)dk-1m 



Re / k '^ + ^ k) d.dk 



kj{uj + n, k) 



dudk 



(12) 

Additionally, we can generalize the empirical stability 
conditions (ESCs) Q to be [H 

f W(dnlmF) + d w ReF - 1< 0, 

\ (d w RcF - l)(d n ImF) - (d w lmF)(d n ReF) > 0, 

(13) 

where dw.n are partial derivatives with respect to W, CI, 
and Re (Im) denote real (imaginary) parts. The ESCs 
(1131) give approximate stability conditions for SSs with 
W > 0. Although empirical, they work well in most 
cases but not all [25|], so should be used carefully. 

Uncoupled distributions. In what follows, we assume 
that the distribution of q, /3, 7 is uncorrclatcd with uj, k: 

G(u, k, q, (3, 7) = g(u, k)h(q, /?, 7), = g(±u + Q, k), 

(14) 

so that the effective lo, k distributions (UJ) become 

I(u, k) = |/|e^ 3 (a;, k), J{w, k) = \J\e^ J g(w, k), (15) 

considerably simplifying all formulae. 
Thus, SCCs dn]) reduce to 



F w (W,SI) = Wcqb4>j/\J\, F n {W,n) 
where 



-Wsin<j>j/\Jl 
(16) 



F W (W,Q) 
F n {W,n) 



dk 
kW 
dk 

kW 

+00 



\k\W 



\k\W 



gl^K 2 R 2 -uj 2 duj, 



{ J ug^duj 



(17) 



[9n ~ 9n 



^lo 2 - k 2 W 2 dcu} 



\k\W 



are the full analogues of -Fr.q in [§]. 

The incoherence stability threshold condition (TT2l be- 
comes 



kg{u + Q, k) 



kg(Q, k)dk 



dudk = - 
2 cos <j)j 



2sin( 



J 



(18) 



while the ESCs (fT3"|) can be simplified using F = 
\J\e l * J (F w +iF n ). Note that, for uncorrelated g(u, k) = 
g(u))v(k), incoherence stability is fully determined by 
(fc) = J kv(k)dk and not by the particular form of v(k), 
as noticed previously for <& = 1, j8< = 7» = (H, [26| . 
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Furthermore, from ([5]), ([6]) and © it follows that 

R=\I\W/\J\,*-$ = <I>i-<I>j, (19) 

System reduction. It is evident that all macroscopic 
properties of the SSs are completely characterized by 
g(w,k), \ J\ and <j)j, while |/|,</>/ serve merely to specify 
Z dHJ). Thus the SSs of systems with different h(q, (3, 7) 
but the same | J|, tfij have same parameters except that Z 
is rescaled and phase-shifted. Therefore, one can obtain 
the parameters of the SSs for ([I]) as 

n = n SK , W = \J\Rsk, R= \I\Rsk, (20) 
where SK denotes parameters of a reduced model 
h I 71 N 

9i=Wi-^p^M»i-0i-^j) (21) 
3=1 

with the same distribution of u)i,ki defined by g(tu,k). 
Although rather hard to prove rigorously (apart from in- 
coherence, for which it follows from (|18p). the stability 
of the corresponding states of (|2"Tj) will also be the same, 
as indicated by ESCs (IT3|) and numerical evidence. 

Hence, for any distribution of 5, /3, 7 obeying (|14l) . one 
can reduce consideration of ([1]) to the much simpler, 
equivalent Sakaguchi-Kuramoto model (l2~Tj) . This ele- 
gant and unexpected result is illustrated in Fig. 1. Such 
reduction, however, concerns only SSs (t — > 00), whereas 
full evolutions Z(t),Y(t) cannot be obtained from (|2T)1) ; 
the microscopic structures of the "equivalent" SSs might 
also be different. Note that, from ®, W < \ J\,R< \I\ 
i.e. the distribution of q, P, 7 imposes an upper bound on 
the SSs' mean field strengths, which they cannot exceed 
however strong the coupling is. 

Glassy states. An oscillator glass is a state where the 
distribution of phases 6i(t) is always uniform and indis- 
tinguishable from incoherence but, in contrast to the lat- 
ter, the oscillators adjust their frequencies and move syn- 
chronously. Relative to the natural frame, the uniform 
phase distribution of incoherence is achieved through 
the phases moving asynchronously (dynamic disorder), 
whereas in the glassy state the phases are frozen in ran- 
dom positions (static disorder). Glassy states are char- 
acterized by vanishing R and nonvanishing W: 

R = 0, W > \I\ = 0, I J| > 0. (22) 

Representing h(q, /3, 7) = hi(P)h 2 (q,j\P), with a (3- 
conditional distribution h 2 (q 7 j\P), ((2"2"|) takes form 

J e if, h 1 (P)d0 = O,f qe^ + ' () h 1 (P)h 2 (q, 1 \P)dqdpd 1 ^ 0. 

(23) 

The first of (|2"3"|) can be satisfied, e.g. if the hi(/3) have 
equal peaks at opposite /? (i.e. separated by tt): 

h i(P) ~ J~] a n [P{P-(f> n , i/ w )+P(/7-<ft w -7r, v n )], ya n ,v n ,(t> n 

n 

(24) 
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FIG. 1. (Color online) Dependences of the SS parameters 
W (green), R (red) and O (black) on coupling K, as cal- 
culated for the original system fTJ (circles) and obtained 
from the reduced system (|21|I . (I20|I (diamonds); solid lines 
show theoretical predictions. The distribution for the origi- 
nal system was: g(tJ, k) ~ S(k — K)[uj 2 + e~" h(q, /?, 7) ~ 
e- (,J - 1)2/0 - 02 [P(/3-7r/8,0.1)+P(/?+7r/2,0.05)]P(7-7r/3,0.2). 
The simulations used iV = 10 5 oscillators and a Runge-Kutta 
6th order method with time-step 0.01 s for 500 s; the values 
are averages over the last 100 s. 



For the second condition, one requires the distribution 
of q, 7 to be specifically correlated with (3. The simplest 
choices are 7 l = -f3 t + <fo (h 2 (q, ~ 6(7 + P - M), 
or qi = cos(Pi + 4>o), but there are many others (e.g. 
h 2 (q, 7 \P)~L(q~ COS p,A q )). 

However, by definition, the glassy state should have a 
uniform marginal phase distribution 

p(6,t) = J f(e,T,t)dT = l/27r. (25) 

As seen from (fTU]) . this is satisfied when h\{P) = l/2ir 
(corresponds to $ty with v n — > 00). Thus, only for uni- 
form p and suitably chosen h 2 {q 1 j\P) are true glassy 
states possible; states satisfying (|22|) but not (|25|) will 
be called pseudoglassy. Note, that the glassy and pseu- 
doglassy states are essentially SSs with particular micro- 
scopic structures, and that they can have fl ^ 0. 

As a simplified picture, the distinctions between the 
states can be understood in terms of a group of people 
doing cyclical exercises, each with their own tempo and 
other parameters. The incoherent state is when every- 
one proceeds independently; synchronized is when they 
all move synchronously, at any given time having identi- 
cal poses; pseudoglassy is when they all exercise at the 
same tempo but remain pairwise in opposite poses; and 
glassy is when they synchronize their tempos, but always 
remain in the independent random poses, thus represent- 
ing a kind of synchronous disorder. Figure [2] illustrates 
these four different states in the natural frame, where 
the population does not move as a whole ((w) =0). In 
the 9, w-plane the only difference between incoherence 
and the glassy state is that, in the latter, phases within 
the glassy cluster (|w| < do not move, remaining 

frozen around random angles P, as seen in 9,P plane. 

Relaxation to an SS. There has been a continuing con- 
troversy about the relaxation of R for glassy states. Some 
authors reported it to be algebraic [10, [l3| , others expo- 
nential [ll|, [12J, [3. Although one should not expect it 
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FIG. 2. (Color online) The different kinds of states, as illustrated by snapshots of positions in the 9, w (upper panel) and 6, /3 
(bottom panel) planes of (the same) 1000 randomly selected oscillators (out of N = 25600). Tiny vertical arrows show the 
movements of the oscillators during 0.25 s. Red dashed lines show boundaries between the clusters and incoherent populations. 
In all cases g(u,k) = L(u, i)5(k-K), h(q,/3,y) = hi(P)6(q-l)6(J3 + j) and (a,b) K = 1, hi(/3) = 1/2tt; (c,d) K = 3, hi(J3) = 
l/2n; (e,f) K = 3, h x [fi) = (l/2)[P(/3 - tt/2, 0.1) + P(/3 + tt/2, 0.1)]; (g,h) K = 3, fti Q9) = P{fi, 0.1). The pictures refer to the 
natural frame ({to) = 0), and all states have fl = 0; for TWs (fi ^ 0) one would observe similar pictures in the rotating frames 
(or add ft to the "speed" of all oscillators in the natural frame). A dynamical version of this figure is available [27l ]. 
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FIG. 3. (Color online) Time-evolution of W(t),R(f) from 
the initial conditions i?(0) = 1 (6i = 0) for the glassy state 
(g(u, k) = S(k - K)L(u, 1), h(q, 13, 7 ) = S(q - 1)6{P + 7 )/2tt) 
for different constant couplings k% — K. The thick gray line 
on the background shows the relaxation of R(t) in the case 
of incoherence (which is independent from K). The inset 
shows the results for a logarithmic ordinate scale. The simu- 
lations used N = 10 6 oscillators and a Runge-Kutta 6th order 
method with time-step 0.01 s. 

to be the same for different KM variants, we find that, 
for ([I]), relaxation is more complex, though close to ex- 
ponential at short and long times, as shown in Fig. [3] 

Surprisingly, for /3-distributions (|24l) the relaxation of 
R to incoherence is fully universal: it does not depend 
on coupling strength, {a„, u n , </>„} in (I24[) or even on the 
form of /i2(S;7|/3). Furthermore, the relaxation to glassy 
and pseudoglassy states occurs in two stages. The phases 
first begin to disorder in the same way as for incoherence, 
so that R decays exponentially with a fixed, coupling- 
independent exponent. At the same time, oscillators are 
also slowly entrained while passing their equilibrium po- 
sitions. When the effective field of the entrained oscil- 
lators, characterized by W, becomes large enough, they 



begin to force unentrained ones to take their positions, 
so the relaxation switches to a faster, coupling-dependent 
exponent and R then decays more rapidly. This "switch" 
occurs sooner for stronger coupling. 

Conclusions. We have generalized the approach of Q 
to a wider case (J]) so that, using (fTTjl . (fT2"j) and (|T3")) . one 
can immediately obtain the macroscopic parameters of 
possible stationary states [28[ . In addition to the partic- 
ular cases included within the framework of [9(, the KM 
([T]) trivially includes the Sakaguchi-Kuramoto model [111, 
as well as more sop histicated cases such as the KM with 
shear diversity [26|, 13CJ, l3l| , allowing one readily to repro- 
duce and extend many previous results [32| . Remarkably, 
the behavior of ([TJ with any distribution h(q,f3,j) (fT4"|) 
in the asymptotic limit t — ¥ oo can be obtained from the 
simple Sakaguchi-Kuramoto model (|2Tj) . (|20|) . 

However, our main finding is that the model ([1} can ex- 
hibit an oscillator glass state, which has very exotic prop- 
erties and which appears when both (p?2")) and ([25)1 are sat- 
isfied. This discovery opens new horizons for KM-related 
investigations, e.g. it will be interesting to see whether 
glassy states can appear when additional complications 
33 39] like time lags are introduced into (JTJ). Finally, 
it should now be possible to create, observe and study 
oscillator glass in real systems, where a variety of novel 
phenomena may be anticipated. For example, if some 
physical quantity can be associated with the weighted 
mean field Y © in laser arrays obeying the KM Q, it 
might be possible to construct a laser with zero intensity 
(R = 0) but other nonvanishing effects (W > 0). Similar 
considerations apply to other KM applications. 

This work was supported by the EPSRC, United King- 
dom (Grant No. EP/100999X1). 
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(|13|) to reproduce and extend many results of [2y| (after 
integration over k,/3 all becomes quite simple). 
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